Possible evidence for two-gap superconductivity in CuxTiSe2 
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We report a comprehensive TF-/iSR study of CuxTiSe2. The magnetic penetration depth was 
found to saturate at low temperature as expected in an s-wave SC. As x is increased we find that the 
superfluid density increases and the size of the superconducting gap, calculated from the temperature 
dependence of the superfluid density, is approaching the BCS value. However, for low values of x, 
the gap is smaller than the weak-coupling BCS prediction suggesting that two superconducting gaps 
are present in the sample. 



The transition-metal dichalcogenides (TMDCs) ex- 
hibit a variety of interesting physical phenomena 
Many of the TMDCs have a charge-denisty wave (CDW) 
modulation in the ground state and many others show su- 
perconductivity (SC) in the ground state. There are few 
examples of systems having both CDW and SC, which 
raises interesting questions regarding the way the Fermi 
surface accommodates these, apparently, competing or- 
ders. The interplay between CDW and SC is re-gaining 
attention, mainly in the context of the newly discovered 
pnictide high temperature superconductors [2|. 

Recently, superconductivity was discovered in Cu in- 
tercalated TiSe2 [3,]. The addition of Cu gradually sup- 
presses the CDW transition temperature and above 4% 
Cu a superconducting phase emerges. Tc increases with 
the Cu concentration up to 4.15K for x=0.08, and de- 
creases for higher concentrations. This compound pro- 
vides an opportunity to study in detail the interplay be- 
tween CDW and SC by simply changing x. 

We report a detailed study of the SC state in CuxTiSe2 , 
using Muon Spin Rotation (/iSR). /iSR allows us to mea- 
sure the magnetic penetration depth. A, as a function of 
the Cu doping, and the temperature dependence of A, 
which is intimately related to the SC gap 

Pressed powder samples of CuxTiSe2 with different Cu 
amounts, ranging from x=0.02 to x=0.081, were prepared 
as described in Ref. Copper concentration was cal- 
culated using X-ray diffraction (XRD) and the known 
c-axis calibration [3]. XRD patterns show no signs of 
foreign phases. 

The resistivity as a function of the temperature for 
all the samples is presented in Fig. [1^. All the sam- 
ples show metallic behavior at high temperatures. The 
x=0.044 sample has, in addition, a broad hump at around 
lOOK. This is related to the formation of the CDW in the 
sample. Above x=0.044 this broad hump cannot be ob- 
served anymore, consistent with the CDW phase ending 
at around 4% Cu. 

A zoom-in on the low temperature part for three of the 
samples is shown in Figure [TId. In Figure [TJ; Tc vs. the 
Cu doping can be seen. The Tc of the x=0.044 sample 
was measured using /iSR, since we cannot measure resis- 
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FIG. 1: (a) The temperature dependence of the resistivity for 
4 samples with various Cu concentration, (b) Resistivity at 
low temperatures, (c) Tc Vs. x. 



tivity below 1.8K. So far, the preparation of CuxTiSe2 
was reported by two groups 0, [sj ; they find very similar 
phase diagrams. In Ref. Q superconductivity emerges 
just above x=0.04. Both groups report x=0.11 to be the 
solubility limit of Cu. We find that the x=0.044 sam- 
ple is superconducting with a transition temperature of 
about 1.5K; samples with lower Cu concentration are not 
superconducting down to 50mK as revealed by /j,SR. The 
maximal Tc is 4.17K, for x«0.08, similar to the maximal 
Tc reported previously 

In Fig. [2] we show thermoelectric-power (TEP) data 
taken at room temperature for the samples. The results 
indicate an increase in the carrier density as Cu is added. 
The Seebeck coeffident is negative for all the samples at 
all temperatures as expected for negative charge carriers 
^lectrons). This is in agreement with previous results [3, 
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FIG. 2: The Cu concentration dependence of the 
thermoelectric-power (TEP) at room temperature. Inset: ZF- 
fiSK asymmetry for the x=0.08 sample at various tempera- 
tures. 



Zero-Field (ZF) /iSR measurements were done at the 
Rutherford Appleton Laboratory, using a dilution refrig- 
erator (DR) and a flow cryostat. ZF /itSR is a very sen- 
sitive local probe of magnetic field. The ZF /j,SR results 
for the x=0.081 sample are shown in the inset of Fig. [2] 
We see a very weak gaussian relaxation for all tempera- 
tures down to 150mK, with no sign of magnetism. This 
suggests that probably each Cu ion contributes one elec- 
tron to the sample, and that the Cu ions are in the Cu+^ 
configuration. Cu++ is magnetic and thus is ruled out. 

Transverse-Field (TF) /iSR measurements were done 
at the Laboratory for Muon Spin Spectroscopy (LMU) 
in the Paul Scherrer Institute (PSI), using the General 
Purpose Surface-Muon (GPS) and the Low Temperature 
Facility (LTF) instruments. 

For the /iSR experiments the samples were mounted on 
silver backing plates. The samples were cooled in a mag- 
netic field of lOOOG from above T^ down to 50mK. The 
magnetic field direction was perpendicular to the muons' 
initial polarization direction. In addition, we measured 
the x=0.081 sample at 1.7K in magnetic field ranging 
from 400G to 1500G. 

In order to analyze the TF-/iSR data we used the fol- 
lowing fitting function: 



A{t) = A,c exp 



1 COS{^^Bintt 



Aas exp(--|^i^) cos{jf,Bextt)- 

(1) 

The first term is the contribution of the sample, while 
the second term is the contribution of the silver holder. 
Asc and Aaq are the initial asymmetries, Bint and B^xt 



FIG. 3: The temperature dependence of the relaxation due 
to the flux-lattice formation for various samples with different 
copper amounts. Inset: asc at the lowest temperature Vs. x. 



are the internal field inside the SC and the external ap- 
plied field, respectively. 7^ = 27r x 135.5342 MHz/T is 
the gyromagnetic ratio of the muon, and 4> is the initial 
phase. There are two contributions to the relaxation of 
the muons in the sample, asc, the origin of which is the 
inhomogeneity of the field due to the formation of the 
vortex-lattice in the sample, and (Tdip, the contribution 
to the dipolar relaxation due to the nuclear moments. 
aAg is the very slow relaxation of the muons in silver. 

The data for each sample over the entire temperature 
range are fitted globally. The initial asymmetries, the ini- 
tial phase, adip and the external field are temperature- 
independent. aAg was measured independently. This 
allows us to reduce the number of free parameters and 
increase the accuracy of the fitting procedure, age is pro- 
portional which, through the London relation, is 
proportional to the superfluid density of the sample 0|. 
For anisotropic materials, with an anisotropy ratio larger 
than five, it was shown that the /iSR line width is pro- 
portional to where A06 is the in-plane penetration 
depth 0- As CuxTiSe2 is very anisotropic [s'], all our 
measurements are sensitive only to the in-plane part of 
the penetration depth. 

The temperature dependence of asc for all the samples 
is presented in Fig. [31 It can be seen that the relaxation 
increases as the temperature decreases, and that it sat- 
urates at low temperatures. In an s-wave SC in which 
there is a fully gapped Fermi surface, a thermally acti- 
vated behavior is expected. This leads to an almost con- 
stant superfluid density at low temperatures. In contrast, 
nodes in the gap lead to a linear temperature dependence 
at low temperatures. Our results are consistent with an 
s-wave order parameter in CuxTiSe2, in agreement with 
previous work [§]. 

In the inset of Fig. |3l we show the low temperature 
value of the muon relaxation versus the Cu content x. 
We find that it increases linearly when Cu is added to 
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the system. In a BCS SC Tc is not related to the super- 
fluid density, in general. Here, we believe that because 
the changes induced by the Cu on the structure of the 
sample are very small, it will have a negligible effect on 
the electron-phonon interaction and on the pairing po- 
tential. Cu increases the amount of free carriers which 
will raise the Fermi-energy and increase the density of 
states at the Fermi energy and will increase Tc. Each Cu 
adds one mobile electron to the system, and that electron 
contributes to the superfluid density below Tc. 

In general, extracting the penetration depth from the 
/iSR line width, a, is not a simple task [l^. For the case 
of a perfect vortex-lattice in a single-gap s-wave SC, there 
is a very good approximation by Brandt based on 
a numerical solution of the Ginzburg-Landau equations. 
This approximation works for k > 5 and for a wide range 
of fields, < B/Bc2 < 1 where k is the Ginzburg- 

Landau parameter: 



(7sc[^ls-^] =4.83-10'' X (1-6) 



1 + 1.21 



A ^ [nm] 



(2) 



where h = B / Bc2 ■ 

A possible self-consitency check is to verify that the 
penetration depth calculated using Eq. [2] is field- 
independent. The penetration depth measured by /iSR is 
calculated by measuring the length-scale associated with 
the decay of the super-currents around the vortex-core. 
This length-scale is expected to be field- independent for 
a single s-wave gap SC [l^. In Fig. we show the 
field dependence of a^c for the x = 0.081 sample, taken 
at T=1.77^. The solid line is a fit to the data using 
Eq. H and Hc2(r = 1.7K) = 0.64T as determined from 
magneto-resitance data, p(H), measured using the same 
sample. We get A = 441 ± 5 nm, and as one can see, 
the agreement is reasonable, but not perfect. A possible 
explanation for this very weak field dependence of A is 
that this sample is not a perfect single gap s-wave SC. 

We use the Gorkov approximation for the tempera- 



ture dependence of Hc2 [I4I to extrapolate IIc2(T) to 



low temperatures. We obtain Hc2(0)=0.85T for the x = 
0.081 sample; we can then estimate the zero-temperature 
coherence-length ^(0) = ^/^^/2^^H^ to be I80A. The 
fact that A does not depend on the field and that we get 
K > 5 indicates that Eq. [2] can be used for analyzing our 
/xSR data, at least for the x — 0.081 sample. Next, we 
used Eq. [2] to calculate A(T) using the data presented in 
Fig. [3^ for three of our samples. It is essential to take 
into account the temperature dependence of Hc2 when 
calculating A. For each sample we measured IIc2 at high 
temperature and used the Gorkov approximation [12!] to 
extrapolate IIc2 to lower temperatures. We found that 
Hc2 depends rather weakly on the doping. The pene- 
tration depth at base temperature for all the samples is 
shown in Fig. 2)3. 
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FIG. 4: (a) The dependance of the relaxation on the external 
magnetic field for the x=0.081 sample, (b) A Vs. the copper 
ratio (x). 



The temperature dependence of A, as already men- 
tioned, is set by the SC gap size and momentum de- 
pendence. A(r)~^ for all the samples saturates at low 
temperature indicating the s-wave nature of the gap. 
The gap size can be extracted from the data using this 
relation tlSl: 



A-' (T) 
A-2 (0) 



= 1-2 



9/ 
dE 



E 



A(T) 



£2 - (A {T)y 



-.dE (3) 



where / is the Fermi function, A(T) is the temperature- 
dependent superconducting gap and the integration is 
over energy measured from the chemical potential. We 
used A(T) = A(0)tanh|l.82 [l.018(^) - l]°'^^|to 

describe the temperature dependence of the gap [14|. 
The data and the fitting curves are shown in Fig. [5] 
The black solids line represent the best fit, and the 
red dashed lines represent the expected temperature 
dependence for the gap size predicted by the weak- 
coupling BCS theory, A(0)=1.764kf,Tc. The gap sizes 
we find are: A(0)=0.60(2) meV, A(0)=0.48(2) meV and 
A(0)=0.095(5) meV for the x = 0.081, x = 0.057 and 
X = 0.044, respectively. 

We systematically get gap values that are lower than 
the BCS weak-coupling values, A (0) = 1.764kB. While 
for the X = 0.081 sample, A/ kBTc=1.67 ± 0.09 is in 
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FIG. 5: Temperature dependence of A(T)~^. The black solid 
lines are best fits to the data using Eq. |3] The dashed red 
lines are the BCS prediction. Figure (a) for x=0. 044(1) sam- 
ple, Figure (b) for x=0. 057(1) sample, and Figure (c) for 
x=0.081(l) sample. 



reasonable agreement with the BCS value, for the lower 
Tc samples we get gap-to-Tc ratios of 1.51 ±0.09 and 
0.7 ± 0.2. The small gap values tell us that the sim- 
ple single s-wave gap picture cannot explain our data. 
A(T) is controlled by thermal excitation of quasi-particles 
across the SC gap. In multi-gap systems or in systems 
with an anisotropic gap A(T) is a measure of the small- 
est gap in the system. Similar results, namely small 
A(0)/Tc, have been found in MgBz fld\ and in NbSe2 
|16| , too. MgB2 has two superconducting gaps as was 
seen, for example, using ARPES [17]. NbSe2, has a multi- 
sheet Fermi-surface [18], with two electron-like cylindical 
Fermi-surface derived from the Nb band, centered around 
the r and K points. In addition, there is a holelike Fermi- 
sheet derived from the Se band around the F point. In 
this system, which is structurally similar to TiSe2, two 
SC gaps were clearly identified on the two hole-like Fermi 
sheets 

Our case is different: We find an evolution of the 
A(0)/Tc ratio as the doping level is changed. Based on 
our data, and the available models for analyzing /xSR 
data, we cannot say if in CuxTiSe2 there are two s-wave 



gaps residing on two different Fermi-surface sheets or one 
anisotropic s-wave gap. 

CuxTiSe2 ARPES data clearly show two Fermi 

sheets, one hole-like pocket centered around the F point 
and one electron-like pocket around the L points. The 
electron-pocket was shown to grow with the Cu concen- 
tration, while the hole-pocket seems less sensitive to the 
amount of Cu. Combined with our result that shows 
that the zero temperature superfluid density increases 
with the addition of Cu, this might suggest that there 
is a small SC gap on the hole pocket the contribution of 
which to the overall superfiuid density is decreasing as 
more Cu is added to the system. 

An alternative scenario is one in which there is a single 
anisotropic s-wave gap, probably on the large electron- 
like pocket, which becomes more isotropic as charge is 
added to the system and the CDW vanishes. Recently 
it was suggested, based on ARPES data in NbSe2 (23 |. 
that at the CDW transition, parts of the Fermi surface 
which are nested are gapped out and excluded from par- 
ticipation in superconductivity. It is not clear to us how 
a situation like that would infiuence the temperature de- 
pendence of the penetration depth, but it might lead to 
an anisotropic SC gap 

In summary, we have performed TF-/iSR experiments 
on a set of CuxTiSe2 samples with different Cu concen- 
trations. The saturation of the muon relaxation at low 
temperatures indicates that the Fermi surface in this sys- 
tem is fully gapped, but the temperature dependence of 
the penetration depth cannot be explained by a single 
isotropic s-wave gap. The data indicate that there are 
probably two gaps in this system the contribution to 
which the superfluid density changes with the Cu con- 
centration. 
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